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Abstract 

It is proven that the center of mass (COM or Kohn) oscillation of a many-body system in a 
harmonic trap coincides with the motion of a single particle as long as conserving approximations 
are applied to treat the interactions. The two conditions formulated by Kadanoff and Baym [l| 
are shown to be sufficient to preserve the COM mode. The result equally applies to zero and finite 
temperature, as well as to nonequilibrium situations, and to the linear and nonlinear response 
regimes. 

PACS numbers: 05.30.-d, 05.60.Gg, 03.75.Kk 
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Interacting many-body systems in confinement potentials are of growinginterest in many 



fields, including ions in Penning or Paul traps 



dusty plasmas 



electrons and 



excitons in quantum wells and dots, e.g. [g], |7j] or ultracold atoms and molecules forming 
Bose condensates, e.g. . In all these systems the problem of collective modes is of high 
interest both, from the fundamental point of view of the iV-particle spectrum and for the 
computation of optical and transport properties, e.g. [g], Q. The normal mode spectrum 
has turned out to be an irnportant and sensitive test of the quality of theoretical methods 



in many-body theory, e.g. [8|, [ill, [12] and of their numerical accuracy. 

Among all collective modes, the center of mass (COM) or dipole mode (Kohn or sloshing 
mode) plays a central role. The dynamics of this mode in an iV-particle system can be 
efficiently analyzed using the familiar operators 

N 

= J2 mVlXk =F iPk, (1) 

k=l 

where Xk and pk denote coordinate and momentum of particle k and, to simplify the notation, 
we restrict ourselves to the one- dimensional case (the extension to 3D is straightforward and 
will be performed below). Application of to the ground state creates an excited state 
with excess energy hVL corresponding to an oscillation of the whole system in the confinement 
potential with the trap frequency f2. In the case of non-interacting particles all oscillate with 
the same frequency, and no deformation of the density profile occurs. Mathematically, this 
statement follows from the Heisenberg equation of motion for the operators C*^, 

showing that, as a function of time, C*^ perform only oscillations without change of the ampli- 
tude. Here Hi = J2k i^~^^xk + ^(^fc)) single-particle contribution to the hamiltonian 
involving the external harmonic potential V{x) = 

It is well known [10] that this behavior is not altered by an arbitrary pair interaction 
w{xi — Xj) between the particles leading to the interacting hamiltonian H = Hi + V12, with 
i<j w{xi — Xj). One readily verifies that [V12, C^] — which is a direct consequence 
of translational invariance of the interaction. Thus, the equation of motion ([2]) remains true 
also for an interacting A^-body system if the interactions are treated exactly which will be 
called "Kohn theorem" . Note that this is not sufficient for invariance of the density profile 
n(r, t) which requires additional restrictions on the initial profile n(r, t = 0) which will be 
discussed below. 
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Obviously, in many situations of practical interest, an approximate treatment of inter- 
action contributions cannot be avoided. The simplest model is to approximate V12 by an 
effective mean field (Hartree approximation), V12 Uh{x) = J dx'w{x — x')n{x'), where 
n{x) is the density. Apparently, Uh does not possess translational invariance and it has 
been stated [Sj] that [Uh,C^] 7^ leading to a modification of the COM mode. This state- 
ment is wrong as can be verified by a lengthy but straightforward calculation . In recent 
years, many improved approximations have been discussed in a variety of fields, including 
the Hartree-Fock approximation, Hartree-Fock-Bogolyubov generalized random phase ap- 
proximation (RPA) [8] or the time-dependent local density approximation [l^, and rather 
sophisticated proofs of preservation of the Kohn mode in linear response have been presented, 



e.g. 



Here, we show that preservation of the Kohn mode can be proven for any approximation 
which conserves particle number and the total momentum. We further demonstrate that 
the result holds independently of whether the system was initially in the ground state, in 
equilibrium or nonequilibrium. Finally, the result equally holds for linear response and for 
the case of arbitrary strong excitation. 

A very compact proof can be obtained using the formalism of nonequilibrium Green's 
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19| . To this end we use the second quantization representation of the 



functions, e.g. 

operators in terms of (fermionic or bosonic) field operators 



c±(X,T) = {^^{x,,x[)^\x[,t[)^{xuh)}^^^, , (3) 



with 7^(xi,x'^) = m Q ^^+^1 ip fi—L-Ji±_ After action of 7 the result is taken at equal 
arguments 1 = {xi,ti) = 1' = [x[,t[) which equal the macroscopic (center of mass) coordi- 
nate and time, X = {xi + x[)/2, T = (ti + t[)/2. Below, the central quantity of interest, 
which takes over the role of in Eq. ([2]), is the ensemble average of ([3]) which is ex- 
pressed by the nonequilibrium Green's function, ihG{l, V) = (T^\l/(l)\l/'^(lO) defined on the 
Schwinger/Keldysh contour C, and denotes time ordering on C [l5|: 



(c±(X,T)) = {7±(xi,XxO^^G(l,l')}^,^^^ , (4) 

where 1"*" denotes the limit from above for the time argument, tf = ti + 0. To simplify 
the notation in the derivations below which apply to fermions and bosons it is helpful to 
introduce G which equals G for bosons and —G for fermions. 



To analyze the behavior of the COM mode of the many-body system as a whole, we need 
to compute the time dependence of the space integrated expectation value 

c±(T)= j dX{c\X,T)). (5) 

This time dependence follows directly from the equations of motion of the Green's functions 
and the adjoint equation (first equation of the Martin Schwinger (MS) hierarchy {igI ) 

{indt,-H^{l)}G{l,l') = inf rf2W^(l-2)Gi2(l,2;l'2+)±5c(l-l'), (6) 

Jc 

{-ih\,~H^{l')]G{l,l') = ihj d2W^(l'-2)Gi2(l,2;l',2+)±5c(l-l'), (7) 

where Hi is the single particle hamiltonian defined above and W{1 — 2) = w{xi — X2) 8c{ti — 
^2), with 5c being the delta function defined on C. 

Note that details of the time contour C and of the matrix structure of G are not important 
for the derivations below. We only recall the physical content of Eqs. (O [7]): while the left 
hand sides describe the single particle dynamics, interaction effects (mean field plus corre- 
lations) are contained in the right hand sides. Here, the central quantity is the two-particle 
Green's function G12 which is a generalization of the standard pair distribution function or 
of the binary correlation operators of the BBGKY hierarchy (the BBGKY hierarchy follows 
from the MS hierarchy in the equal time limit, e.g. [3]). Practically all many-body approx- 
imations can be formulated in terms of Gi2- For example, the Hartree Fock approximation 
follows by substituting Gi2(l, 2; 1', 2') ^ ^(1, 10^(1, 1') ± ^(1, 2')G'(2, 1'). Finally, the ad- 
vantage of the formalism of Eqs. ([6l [7]) is that it allows to generalize all approximations 
known from ground state many-body theory, including diagram expansions, to arbitrary 
nonequilibrium situations. 

We now proceed to compute the time dependence of c^(T), Eq. ([5]), and show that one 
exactly recovers equation ([2]). In particular, we show that the contributions from 6*12, i.e. 
from the interactions, to the dynamics of c^{T) vanish exactly if (A): the Green's function G 
obeys the two mutually adjoint equations ([6l[7]) and (B): for G12 a conserving approximation 
is used which only requires Gi2(l, 2; 1', 2') = (712(2, 1; 2', 1'). These are exactly the two 
criteria of Baym and Kadanoff known to be sufficient for conservation of density (continuity 
equation), momentum and total energy 0,0. 

The time derivative of c^(T) is computed directly using the definitions (jlJlS]) and taking 
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the difference of Eqs. dH [7]) 

d 



dT 

- dX \ ^^{xi,xi 



dX {7±(xi, xv) {dt, + \,)inG{l, 1')}^, 



=1+ 



h (5x1 - 



2m 



- [V{x{) - V{x[)] 



G(l,l') 



+ i/;.y dX dx2 |7^( 



1' = 1H 



Xi, Xi/j — X2) — W[Xi' — X2) 



(8) 



where the terms 6c cancel and, in the interaction term, the time integration has been carried 
out. We now apply the operator 7^ under the integrals and obtain the single particle 
contribution dcf/dT to (jH]): 
d 



dT 



ctiT) 



dX 



d - d 
mVtX ^ h — 



X 



X 



(9 - (9 



- / dX 



± 



2m 

V\xi) + V\x'^ 



[V{x,) - V{x[)] 
d - d 



i'=i+ 



i'=i+ 



(9) 
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where a partial integration over X has been performed taking into account that G vanishes 
for |X| — > 00. For a harmonic confinement potential V follows dcf/dT = ±zfic^. 
Consider now the interaction contribution to (Mi: 
d 



^cf2(T) = J dXdx2^ 



m^7X =F h — 



=1+ 



X [w{xi - X2) - w{xr - X2)] 6-12(1, 2; 1', 2+) 1 

J 1'= 

Contributions proportional to VlX vanish because the potentials cancel for The 
remaining contribution involves derivatives of the interaction potential, and the integrand 
can be transformed according to 

{ [d^,w{xi - X2) + d^^,w{xr - X2)] G^il, 2; l', 2+)}^,^_^^ 
= {2w'(X -0:2)^12(1, 2; 1', 2+)}^,^^^ . 

Since the force —w'{x) is an odd function of the argument, the expression in parantheses 
vanishes under the integral over the coordinates of both particles if Gu is even. But this 
is exactly condition (B). Thus, symmetry of Gu with respect to the arguments of particle 
1 and 2 is sufficient for vanishing of the time derivative dcf2/dT. Summarizing the above 
results for ci and C12, we obtain for the sum ([8]) 

d 



dT 



(10) 



Thus, we exactly recover the dynamics of Eq. ([2]) with the only assumption that conditions 
(A) and (B) are fulfilled. Note that our derivation did not involve information on the type of 
pair interaction w and is thus valid for arbitrary pair potentials. Also, no specific ensemble 
had to be specified for the average (in the Green's function), therefore, the derivation is valid 
for A^-body systems originally in the ground state, at finite temperature or in nonequilibrium. 

As was shown by Baym [isl conditions (A) and (B) are equivalent to existence of a 
functional $ such that 1') = 6^[G]/6G{1' ,1). This means we may conclude that any 
derivable approximation" for the selfenergy fulfills the Kohn theorem. This result has 
important implications not only for Green functions theory but also for other classes of ap- 
proximations, such as the ones used in time-dependent density functional theory (TDDFT). 
It was recently proven j20] that any exchange correlation functional of TDDFT which is de- 
rived from a $— derivable approximation satisfies the so-called Zero-Force Theorem. Thus, 
conditions (A) and (B) are suitable to construct TDDFT approximations which preserve 
the Kohn mode [l^ . 

A number of generalizations are evident. Our result is trivially generalized to the 3D 
(2D) case with a confinement of the form Z^i^i^^f^jf, where M = 3(2). In this case, there 
exist three (two) independent COM oscillations, and the derivations can be repeated sep- 
arately for each of the three (two) expectation values q(T) which now involve Qi instead 
of Q. Also, while our derivation was performed for a single component system, general- 
izations to multicomponent systems, including Bose condensates, superfiuid Fermi gases or 
mixtures of fermions and bosons, should be straightforward. Furthermore, the derivations 
presented in this work are readily generalized to systems in an additional one-particle po- 
tential U{r,t) = Uoit) + Ui(t)r, such as an electromagnetic dipole potential (f/o and Ui are 
spatially homogeneous), since this only leads to a replacement ^(r) — > V{r) + f/(r, t) in the 
single particle hamiltonian Hi. In that case, the equations of motion ([2]) of the operators 
are changed to ^C*^ = ±zfiC'='= ± iNUiif). From this it follows that the center of mass 
R = j;jT.kXk = {C^ + C-)/{2niNn) obeys the equation R + fi^i? = -Ui{t)/m, which is 
just the classcial equation of the forced harmonic oscillator. Again, this result is not altered 
by interactions as long as they are treated within approximations obeying conditions (A) 
and (B). 

Our result can also be directly extended to other types of collective modes. As an example, 
we mention the cyclotron oscillations of charged particles in a homogeneous magnetic field 
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(0,0,5). Then, the operators (7=*= describing the modes have to be replaced by |lO| P± = 
Px ± iPy with P = {px,Py + ^,Pz)- Again, it is known that the cyclotron rotation with 
ujc = eB/{mc) (which is excited by action of on the ground state) is not altered by 
pair interactions if they are treated exactly jlo|. The question of invariance of the cyclotron 
rotation under an approximate treatment of the interaction can be analyzed exactly as above. 
Instead of the operators d^^ — dx'_^ now there will be linear combinations of different spatial 
components of the momentum operator which again lead to vanishing of the interaction 
contributions dcu/dT, for any conserving approximation. 

Finally, let us come back to the question under what conditions the density profile n(r, t) 
will remain unchanged (in the frame oscillating with the center of mass). It is easily verified 
that invariance of n(r, t) in an interacting system is a direct consequence of the fulfillment 
of Eq. ([2]), if initially the system was prepared in an eigenstate of or in a (correlated) 
thermodynamic equilibrium state 

In summary, we have shown that the center of mass oscillation (Kohn or sloshing mode) 
of interacting particles in a harmonic potential remains unchanged also in an approximate 
treatment of the many-body problem, for any conserving approximation. This result is of 
fundamental importance for many fields of physics where trapping potentials are being used, 
including ultracold ions, dusty plasmas, nanostructures or quantum gases. For theoretical 
treatments, our result allows to single out the scope of physically reasonable approximations, 
which are just the $— derivable approximations of Baym. At the same time, it predicts that 
the shape of a density distribution of interacting classical or quantum particles in a harmonic 
trap should remain unchanged in time, if initially only the center of mass oscillation is excited 
and the initial state was either an eigenstate of the hamiltonian or equilibrium. 
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